**FULL TITLE** 

ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** 
**NAMES OF EDITORS** 



Missing Halo Baryons and Galactic Outflows 



Romeel Dave 
University of Arizona 

Abstract. We present predictions for galactic halo baryon fractions from 
cosmological hydrodynamic simulations with a well-constrained model for galac- 
tic outflows. Without outflows, halos contain roughly the cosmic fraction of 
baryons, slightly lowered at high masses owing to pressure support from hot 
gas. The star formation efficiency is large and increases monotonically to low 
masses, in disagreement with data. With outflows, the baryon fraction is in- 
creasingly suppressed in halos to lower masses. A Milky Way-sized halo at z = 
has about 60% of the cosmic fraction of baryons, so "missing" halo baryons have 
largely been evacuated, rather than existing in some hidden form. Large halos 
(;> 10 13 M Q ) contain 85% of their cosmic share of baryons, which explains the 
mild missing baryon problem seen in clusters. By comparing results at z — 3 and 
z = 0, we show that most of the baryon removal occurs at early epochs in larger 
halos, while smaller halos lose baryons more recently. Star formation efficiency is 
maximized in halos of ~ 10 13 Mq, dropping significantly to lower masses, which 
helps reconcile the sub-L» slope of the observed stellar and halo mass functions. 
These trends are predominantly driven by differential wind recycling, namely, 
that wind material takes longer to return to low-mass galaxies than high-mass 
galaxies. The hot gas content of halos is mostly unaffected by outflows, showing 
that outflows tend to blow holes and escape rather than deposit their energy 
into halo gas. 



1. Introduction 

Everywhere we look, baryons are missin g. Global l y, ob servations only account 
for around half the baryons today (e.g. iFukugital 20041 ). Simulations indicate 
that the missing baryons are contained in intergalactic g as at 10 5 < T < 1 7 K, 
called the Warm-Hot Intergalactic Medium (WHIM; e.g. iDave et al.l[200ll ). Ini- 
tially, O vi absorbers seen i n quasar spectra were thought to be collisionally- 
ionized tracers of WHIM gas ( Tripp. Savage. &: Jenkins1l2000l ). But recent simu- 
lations and observations suggest that most intergalactic O vi is actually photo- 
ionized ( Qppenheimer &: Dave 2008bl . and refs ther ein), and so the WH IM must 
be traced using higher ionization lines (e.g. O vil; iNicastro et aLll2005l ) . 

Baryons are also missing on galactic scales. Current hierarchical structure 
formation models predict that baryons do not substantially decouple from dark 
matter until well inside of halos, hence the expectation is that halos should 
contain roughly the cosmic fraction of baryons. However, dynamical model- 
ing of the Milky Way's disk and halo reveals that it contains at best half 

of its cosmic share of baryons in stars and cold gas ( Dehnen Sz Binnevl 19981; 

So mmer-Larsen fe Dolgov 200ll ); the same is true of M31 ( Klvpin. Zhao. Somerville 
20021 ). This is the "Missing Halo Baryons" problem. Either a substantial portion 
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of the ~ galaxies' halo baryons are in some heretofore hidden form, or else 
there has been a sizeable exodus of baryons into the intergalactic medium (IGM). 
Many models have advocated the former possibility, suggesting by analogy with 
the WHIM that the missing halo gas is in s ome warm-hot component (e.g. 
Fukudta &: Peeblesl l2006 ; Sommer-Larsenll2006l ) , or cold clouds (jMaller Sz Bullockl 



20041 b However, there is lit tle observational ev i dence for m assive coronae of hot 



gas around typical spirals ( Benson et al.l l200(il ; IWangl 120071 ) . and furthermore if 
all halos containe d such coronae, this would substantially overpredict the soft X- 
ray background ( Penlll999l ; IWu. Fabian. fe Nulsenll200ll ). Meanwhile, the mass 
in cold clumps as tra ced by high velocity clouds is unlikely to be sizeable ev en in 
optim istic scenarios ( Blitz et al. 19991 ; Sommer-Larsenl 120061 ) . In contrast, Silk 
(2003) developed an analytic model that argued for a substantial fraction of 
baryons being removed by galactic outflows. In fact, Silk presciently predicted 
that such outflows could solve a surprisingly wide range of current dilemmas in 
galaxy formation, which our simulations have largely confirmed. 

Clusters also seem to have a missin g baryon problem. In a recent census by 
IGonzalez. Zaritskv. Zabludoff I ( 20071 ). baryons make up sa 13% of t he mass 
of the cluster, while the latest WMAP-5 results ( Hinshaw et al.l 120081 ) favor a 
cosmic mean value of ~ 17%. This is a small but persistent discrepancy which 
could have implications for using clusters as probes of precision cosmology. 

In these proceedings we examine the baryonic content of galaxy halos in cos- 
mological hydrodynamic simulations with and without galactic outflows. The 
outflow model implemented in our simulations is unique in that it matches de- 
tailed properties of a wide rang e of outflow-related observables, inc l uding IGM 
enrichment at various epochs (iQppenheimer fc Pavel 120061 . l2008bl . l2009t >. the 
galaxy mass- metallicity relation (Finlator &: Dave 20081 ). early galaxy luminos- 
ity functions ( Dave. Finlator. &: Qppenheimerll2006h . and int ragroup gas enrich- 
ment and entropy (|Dave. Oppenheimer. &: Sivanandamll2008l ) . Even though our 
modeling of outflows is parameterized and heuristic, these successes suggest that 
it plausibly moves mass, metals, and energy on large scales in a manner consis- 
tent with the real Universe. 



2. Halo Baryon Fractions 



We run sim ulations using our mod i fied ve rsion of Gadget- 2 ( Springell2005l ). as 



described in IQppenheimer fc Pavel (|2008al b The runs here employ 2 x 256 3 par- 



ticles in volumes of 32 h 1 Mpc and 64/t 1 Mpc, with our WMAP3-concordant 
d-series cosmology (see IQppenheimer & Pavel l2008blb We identify halos using 
a spherical overdensity algorithm (see lKeres et al.ll2005l ). and consider only "re- 
solved" halos with masses Mh > I28(m d + m g ), where m d = (1.6, 12.7) x 1O 8 M 
and m g = (0.34, 2.72) x 1O 8 M0 are the dark matter and gas particles masses 
for the (32, 64)/i _1 Mpc volumes, respectively. We run two versions of these 
simulations: One with no winds, and one with our favored momentum-driven 
(M-P) wind model that matches a range of data as mentioned above. Outflows 
are implemented in a probabalistic Monte Carlo fash ion by giving kicks to gas 
particles; details are in IQppenheimer fc Pavel (|2008al b 



The fraction of baryons within halos are shown in the top panel of Figure [TJ 
In the no- wind case (blue points, upper swath), halos of Milky Way's size and 
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Figure 1. Top: Baryon fraction (relative to cosmic mean) as a function of 
halo mass, for simulations at z = (left panels) and z — 3 (right), with- 
out winds (blue, upper points) and with momentum-driven winds (green, 
lower points). Dotted line shows the cosmic baryon fraction assumed in our 
runs. Middle: Stellar fractions. Bottom: Hot gas (T > 10 4 ' 5 K) fractions. 
In each panel, the red line shows a running median with la dispersion, and 
the two groups of points for each model show results from our 32ft, _1 Mpc 
and 64/i -1 Mpc volumes. Outflows drive out substantial amounts of baryons 
particularly from low-mass systems, which helps regulate star formation to 
observed level. In large systems most of the mass loss has already occured by 
z = 3, while in small systems the mass loss occurs later owing to differential 
wind recycling. Hot gas fractions are a strong function of halo mass, but are 
not sensitive to outflows. A Milky Way-sized halo today contains about 60% 
its cosmic share of baryons, of which half is in the form of hot gas, one-third 
in cool gas, and one-sixth in stars. 
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below have, as expected, roughly their cosmic share of baryons (shown by the 
dotted line at fi&/fi m = 0.044/0.25 = 0.176). In detail, adiabatic contraction 
causes the baryon fraction to typically be slightly greater than the cosmic mean. 
This shows that without any strong feedback, baryons do indeed mostly trace 
dark matter on halo scales in <^ L* galaxies. 

In m ore massive halos (M^ ^ 1O 13 M ), the increasing predominance of hot 
halo gas ( Keres et al.ll2008l . see also bottom panels of Figure [Q) causes pressure 
support that pushes baryons farther out than the dark matter, and the baryon 
content is lowered. At the highest masses probed by z = (M h ~ 1O 14 M ), 
this reduction reaches 6%, which is not trivial but still insufficient to explain 
observed cluster baryon fractions. 

Now we consider the M-D wind runs, shown by the lower (green) swath 
of points. Outflows have a dramatic impact, increasingly so to smaller halo 
masses. Already by z = 3, typical halos (Mh ~ 10 11_12 Mq) have had their 
baryon fraction reduced by 20 — 30% compared to the cosmic mean. At z = 3 
there is only a mild trend with halo mass, but by z = the trend becomes much 
stronger: Milky Way-sized halos (Mh ~ 10 12 M Q ) have lost roughly 40% of their 
baryons, while the smallest halos we probe at Mh ~ 1O 1O ' 5 M0 have lost 60%. 
This indicates that missing halo baryons have mostly been ejected by outflows. 



Outflows have a noticeable impact on large halos as well. iDave. Oppenheimer. Sz Sivanandam 
( 20081 ) showed that outflows add substantial entropy on poor group scales, and 
this translates into increased pressure support over the no-wind case. With out- 
flows, the most massive (Mh ^ 1O 13 M0) halos now contain only 85% of their 
cosmic share of baryons, and the trend is essentially flat with halo mass. This 
value is in better agreement with obser vational estimates ( Vikhlinin et al. 20061 : 
Gonzalez, Zaritskv. Zabludoff ll2007l ). 



The falling baryon fraction in low-mass halos is a consequence of wind 
recycling in our outflow model, i.e. the re- accretion of previously ejected wind 
material. In Oppenheimer Pavel ( 2008al ). we found that small galaxies push 
their winds out farther relative to large galaxies (owing primarily to the fact that 
they reside in less dense surroundings), and as a result the time for winds to be 
re-accreted onto small systems is longer. At early epochs, recycling is relatively 
unimportant, because there hasn't been sufficient time to permit substantial 
re-accretion. The weak trend seen at z = 3 owes to the fact that our mass 
loss rate scales inversely with velocity dispersion in our momentum-driven wind 
scalings. By z = 0, however, the longer recycling times in smaller systems 
results in much more cumulative mass loss from small halos relative to larger 
halos. Hence differential wind recycling is the dominant driver in establishing 
the baryon fraction trends with halo mass. As we will see next, it is also critical 
for regulating the star formation efficiency in smaller systems. 



3. Stellar Baryon Fractions 

We now separate baryons within each halo into three phases: Stars, cool gas 
(T < 10 4 ' 5 K), and hot gas (T > 10 4,5 K). Star- forming gas is included in cool 
gas. The middle and bottom panels of Figure Q] show the baryon fractions in 
stars and hot gas, respectively; the remainder is in cool gas (not shown). In this 
section we examine stellar baryon fractions. 
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With no outflows, there is a serious overcooling problem: Star formation 
is far too efficient. In sub-L* halos, the fraction of baryons in stars approaches 
90%, while observed small galaxies have stellar fractions under 10%. Moreover, 
the trend is wrong; with no winds, smaller galaxies have a higher baryon fraction 
in stars, while the shallow slope of the stellar mass function below relative to 
the halo mass function indicates that the stellar mass fraction should be lower 
in small halos. At group scales, the stellar baryon fraction is down to 25% of 
th e cosmic mean, but that is still too high compared to data as pointed out 
in lDave. Qppenheimer. Sivanandaml ([2008). Of course, these results are not 



surprising, as it is well known that outflows, particularly at early epochs, are 
required to suppress overcooling (e.g. Dave. Finlator. Qppenheimer 20061 ). 

With outflows, the situation improves dramatically. Stellar fractions are 
reduced to a maximum of ~ 20% for ~ 10 13 Mq halos, falling to either smaller 
or larger halos. A Milky Way-sized halo now has a 12% stellar fraction (relative 
to fif,), which is in general a greement with observations of comparable disk 
galaxies ( Hammer et al.l 120071 ) . The stellar fraction drops towards lower halo 
masses, suggesting a flatter stellar mass function compared to the halo mass 
function, as observed. The slope of the z = fit (red line) is approximately 
dlog f*/dlog Mh « 0.5 for Mh < 5 x 10 12 M Q , which when combined with the 
halo mass function slope of rj —2, yields a stellar mass function slope of ~ —1.5. 
This is in general agreement with observatio ns of the sub-L* stellar mass function 
slope (jBaldrv. Glazebrook. &: Driver! l2008l ). though still too steep; in fact, the 
faint-end slope turns out to be even shallower, as will be shown in a forthcoming 
paper. Note that the stellar fraction slope at z = 3 is less steep, reflecting the 
shallower slope in the overall halo baryon fraction (top right panel). Our models 
naturally yield a flattening of the faint end slope with time. 

The qualitative trend of the star formation efficiency having a maximum at 
some mass and dropping fairly rapidly to low masses, along with the faint end 
slope evolution, is a direct consequence of differential wind recycling. Without 
it, the star formation efficiency would continue to increase to small masses, as 
in the no-wind case. Hence models of galaxy formation must not only include 
outflows, but must also track the dynamics of wind material on its journey 
through the IGM to properly capture its impact on galaxy evolution. 



4. Hot Gas Baryon Fractions 

Turning to the hot gas (T > 30, 000 K) baryon fractions (bottom panels of 
Figure [1]), there is a strong trend of hot gas fraction increasing with halo mass, 
which is mostly independent of redshift and whether or not outflows are included. 
In the no-wind case, the hot gas fraction exceeds the stellar fraction (which 
comprises the vast majority of condensed baryons) at z = for > 10 12,5 M^. 



This transition occurs at a somewhat larger mass than found in iKeres et al.l 
(120081 1. likely because our runs include metal-line cooling. 

With outflows, the fraction of halo gas in hot form is more substantial, 
exceeding the stellar fraction in the outflow case at all masses probed here. 
This is due to the suppression of star formation, rather than an increase in 
the amount of hot gas. Relative to the baryon fraction, the hot gas fraction is 
about one-third of all halo baryons at the smallest masses probed, and increases 
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to three-fourths or more for ^ 1O 13 M0. This indicates that there is still 
a substantial reservoir of baryons in a difficult-to-detect phase within typical 
galaxies, and indeed X- ray observations do i ndicate a warm-hot corona around 
a nearby spiral galaxy ( Pedersen et al. 20061 ). However, the amount of hot gas 
we predict is still a factor of two to three smaller than in models that place the 
majority of missing halo baryons into this phase. This roughly translates into up 
to an order of magnitu de reduction in the predicted X-ray flux (as compared to 
e.g. iBenson et al. 20001 ) . Hence the typical non-detection of X-ray halos around 
spiral galaxies should not be surprising, but deeper observations with future 
X-ray telescopes should uncover this phase more ubiquitously. 

It is noteworthy that our outflows do not significantly increase the amount 
of hot gas into halos (as seen by comparing to the no-wind case). Winds in our 
models tend to escape halos without depositing a significant amount of energy 
along the way via shocks. This is cont rary to typical assu mptions in analytic 
or semi-analytic models of outflows (e.g. iDekel Sz Silklll986l ). and simply reflects 
the fact that in realistic three-dimensional models outflows prefer to blow holes 
rather than share their energy with ambient gas. This is also consistent with 
X-ray observations of hot gas around galaxies that indicate the y are radiating a 
small fraction of the supernova energy input (e.g. Wand [20071 ) . Another factor 
is that outflows are typically quite enriched, so that the gas cooling times are 
short, and hence even if outflow material is heated it condenses out quickly. At 
low-z this resul ts in a "halo fountain" that may be the o rigin of compact high 
velocity clouds (jWangl [20071 ; lOppenheimer Pavel l2008al ) . 



5. Summary 

We have examined halo baryon fraction in cosmological hydrodynamic simula- 
tions, comparing models with and without galactic outflows. Our outflow model 
is heuristic, but is well-constrained to match a variety of galaxy and IGM obser- 
vations at a range of epochs. We find that these same outflows may help resolve 
some puzzles regarding the baryonic content of galactic halos. 

For low mass halos, outflows remove a significant portion of halo baryons 
by z = 0. The fraction of ejected baryons increases sharply to lower masses, 
so that 1O 12 M halos have lost 40% of their baryons, and lO ia5 M halos 60%. 
The trend is not as pronounced at z = 3, showing that the longer wind recycling 
times in small galaxies (i.e. differential wind recycling) plays a critical role in 
suppressing the baryon fractions in small halos between z = 3 — > 0. This may 
explain why the faint-end slope of the luminosity function becomes shallower 
with time. By today, this produces a peak in star formation efficiency at ~ 
10 12 ' 5_13 Mq, above and below which the efficiency falls. Meanwhile, hot gas 
fractions are mostly unaffected by outflows, showing that outflows do not deposit 
much of their energy into halo gas, instead preferring to blow holes and escape. 
A Milky Way-sized halo today has about 60% of its cosmic share of baryons, 
half of which are in hot gas, one-third in cool gas, and the remainder in stars. 

For high mass halos (M/j ^ 1O 13 M0), there is a mild suppression of baryon 
fractions owing to larger pressure support from an increasingly substantial hot 
gaseous halo. Without winds, the suppression is fairly small (~ 5%). In our 
outflow run, conversely, the suppression is much more substantial, ~ 15%, inde- 
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pendent of mass, which agrees better with data. This demonstrates that outflows 
impact even massive halos today. The key is that most of the ejection occurs at 
early epochs, when those halos were much smaller. 

Our implementation of o utflows is primarily constrain ed to match IGM en- 
richment at z ~ 2 — 4 (e.g. lOppenheimer &: Dave 2006). The fact that this 



model naturally yields observationally-consistent results for halo baryon frac- 
tions at z ~ is highly encouraging, and represents another significant success 
for our outflow model based on momentum-driven wind scalings. It may be 
that we are approaching a heuristic understanding of how outflows operate on 
extragalactic scales, even though the detailed mechanisms by which outflows are 
driven out of galaxies remain unclear. 

Acknowledgments. RD thanks Shardha Jogee and the organizing com- 
mittee for an excellent meeting. The simulations were run on UA's SGI cluster. 
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